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^— N ■ The Stephani universes that can be interpreted as an ideal gas 

^+ I evolving in local thermal equilibrium are determined. Five classes of 

^^ ■ thermodynamic schemes are admissible, which give rise to five classes 

fi I of regular models and three classes of singular models. No Stephani 

^^i universes exist representing an exact solution to a classical ideal gas 

t , I (one for which the internal energy is proportional to the temperature). 

bi)' But some Stephani universes may approximate a classical ideal gas at 

first order in the temperature: all of them are obtained. Finally, some 
features about the physical behavior of the models are pointed out. 
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1 INTRODUCTION 

The second member of Einstein equations describes the energy content of 
the medium, but not the medium itself. So, exact solutions to Einstein 
equations correspond to prescribed forms of energy tensors, but they do not 
fix the nature of the particular (class of) physical fluid(s) that creates them. 

Generically P, the physical interpretation of an exact solution to Ein- 
stein equations, needs more than that of their energy quantities, i.e. more 
than that of the quantities that appear in the energy tensor: it is also nec- 
essary to 'close the context'. Let us see what this means. 

A non vacuum physical space-time is the history of one of the possible 
evolutions of a physical medium. But there are in general many physical 
media that can have in common one or some particular evolutions (for exam- 
ple, many fluids admit a rigid and static evolution with same energy density 
an pressure). In fact, from the evolution point of view, the definition of a 
medium is tantamount to the knowledge of the set of all its possible physical 
evolutions. Of course, this knowledge does not need to be extensive, it is 
sufficient that it be intensive, i.e. that it be done by the set of equations 
whose space of solutions is precisely the whole set of all its possible physical 
evolutions. 

But according to a universal determinism requirement in Physics, such 
a set of equations has to be (causally) closed, i.e. has to admit a well posed 
initial value problem. 

In general the conservation equations imposed by Einstein equations 
on the energy tensor (energy conservation equations) do not constitute a 
closed system of equations for the energy quantities. This is why we say 
that their physical interpretation is not sufficient to define the medium. 
To 'close the context' means here just to choice a closure for the energy 
conservation equations, i.e. to complete these last ones with some other 
physically meaningful equations in such a way to obtain a closed system. 

The closure of the energy conservation equations usually 01 obliges to in- 
troduce new quantities, related at least in part to the Duhem-Gibbs balance 
equation, and called by this fact thermodynamic quantities. 

Thus, for a perfect energy tensor, 

T = {p + p)u®u-pg , (1) 

where p is the energy density, p the pressure and u the (unit) velocity of 
the fluid, the usual closure .S; corresponds to the hypothesis of local thermal 
equilibrium (also called standard thermodynamic scheme). The local thermal 
equilibrium evolution of a fluid is thus determined by the following four 
conditions: 

• Energy-momentum conservation: V • T = 0. 



• The energy density p is decomposed in terms of the matter density r 
and the specific internal energy e: p = r(l + e). 

• The equation of conservation of matter is required: V • (ru) = 0. 

• The thermodynamic quantities temperature T and specific entropy s 
are related by equations of state compatible with the thermodynamic 
Duhem-Gibbs relation: Tds = de + pd[\/r). 

It is clear that this standard thermodynamic closure is a strict constraint 
on the whole space of solutions to the Einstein Equations for perfect energy 
tensors, i.e. not all the solutions of this last space can be interpreted as 
being local thermal equilibrium evolutions of a perfect fluid |^ . 

Observe that this constraint notably increases the number of variables 
because it involves not only the starting energy quantities (it, p, p) but also 
the new thermodynamic ones (r, e, T, s). So that, an important question 
arises: is it possible to know, by the sole inspection of the energy quantities of 
a perfect energy tensor, if it describes the local thermal equilibrium evolution 
of a perfect fluid? 

As we showed some years ago j^, the answer is positive. Its precise 
formulation (denoting by u{x) the directional derivative, with respect to u, 
of the quantity x) is given by the following 

Theorem x (local thermal equilibrium): A divergence-free perfect en- 
ergy tensor evolves in local thermal equilibrium if, and only if, the space-time 
function x = u{p)/u{p) depends only on the variables p and p : dxf\dpf\dp = 
0. 

This result is strikingly interesting because it shows that, in spite of its 
name, of its historic origin and of its usual conceptualization, the notion of 
local thermal equilibrium for a perfect fluid is a purely hydrodynamic, not 
thermodynamic, notion 0. 

The function x is called the indicatrix of local therm,al equilibrium. To 
every physical fluid in local thermal equilibrium, whatever its evolution be, 
corresponds a sole, particular, indicatrix (particular equation of state) which 
physically represents the square of the sound velocity in this fluid |7|. So 
that the direct application of theorem x iiot only allows 

i) to detect, among the known families of solutions to Einstein equa- 
tions for a perfect energy tensor, those sub-families that describe evolutions 
in local thermal equilibrium of physical fluids, but also 

ii) to separate these last sub-families in sub-sub-families that describe 
different (local thermal equilibrium) evolutions of the same physical fluid. 

Conversely, theorem x generates the first general method to know wether 
a solution for a perfect energy tensor is the evolution of a previously given 
fluid: it is sufficient to verify if the indicatrix function x of the solution 



coincides with the square of the thermodynamic expression of the sound 
velocity of the given fluid expressed as a function of the variables p and 
p. In we pointed out these applications and gave a concrete example 
studying the thermodynamic class II Szekeres-Szafron space-times. 

In the present paper we consider, on one hand, the paradigmatic ideal 
gas as the previously given particular fluid and, on the other hand, the 
Stephani universes jSj as starting family of solutions to Einstein equations 
for a perfect energy tensor, both because of their relative generality and 
their frequent applications in cosmology and stellar interiors. 

The Stephani universes admitting a generic thermodynamic scheme were 
obtained by Bona and Coll [H] and later revisited by Krasiiiski et al. ^01 • In 
a recent paper, Sussman ^J has analyzed a family of spherically symmetric 
Stephani universes that may be interpreted as either a classical mono-atomic 
ideal gas or a matter-radiation mixture. 

Here we are interested in the deductive and systematic determination of 
all Stephani universes that represent the evolution in local thermal equilib- 
rium of a generic ideal gas and, among them, of those that admit a first 
order contact in the temperature with the classical ideal gas |12j . 

The main results of this work are based on previous ones: firstly, on 
our hydrodynamic characterization of an ideal gas in local thermal equilib- 
rium pL3 , that we shorten here in section 2; secondly, on the above cited 
work by Bona and Coll [H] about the thermodynamic schemes in Stephani 
universes, that we summarize in section 3. In section 4 we apply an algo- 
rithm introduced in section 2 to a canonical form for the thermodynamic 
Stephani universes presented in section 3 in order to obtain: i) a theorem 
(theorem 1) that characterizes the generic ideal gas Stephani universes and 
ii) the associated indicatrix function and five classes of compatible ideal gas 
thermodynamic schemes. In section 5 we integrate, up to a Priedmann-like 
equation, the ideal gas equations obtained in theorem 1; this integration 
leads to different models with thermodynamics of a specific class. In section 
6 we study the Stephani ideal gas models that approximate a classical ideal 
gas at first order in the temperature. Finally, section 5 is devoted to analyze 
how complementary physical requirements, like energy and compressibility 
conditions, restrict the domain where the models have a good physical be- 
havior. A part of the results of this paper were stated without proof in 

EI- 

2 LOCAL THERMAL EQUILIBRIUM EVOLUTION OF 
AN IDEAL GAS 

As already mentioned. Theorem x provides an easy and general method (in 
fact the first known one) to answer if whether or not a solution to Einstein 



equations for a perfect energy tensor corresponds with a given thermody- 
namic fluid: to check that, as functions of p and p, the indicatrix function 
of the solution and the square of the sound velocity of the fluid coincide. 

Elsewhere ^3j we have given the hydrodynamic characterization of the 
sound velocity for some representative classes of fluids. We present here in 
two propositions the part of these results about ideal gases that we will use 
in this paper. 

Remember that a perfect energy tensor is said barotropic if p and p are 
not independent, dp A dp = 0, and that it is said isoenergetic if u{p) = 0. As 
we shall see in the next section, for the class of space-times analyzed here, 
the Stephani universes, the barotropic case leads to the already well known 
Friedmann-Robertson- Walker models, meanwhile the isoentropic case is not 
possible because of the non vanishing expansion of these fluids. 

A (generic) ideal gas satisfies the equation of state p = krT. For it we 
have shown the following hydrodynamic characterization |13) : 

Proposition 1 The necessary and sufficient condition for a non barotropic 
and non isoenergetic divergence-free perfect energy tensor (^ to represent 
the local thermal equilibrium evolution of an ideal gas is that its indicatrix 
function x = u{p)/u{p) be a non identity function of the variable n = p/p: 

dxAdTT = 0, X ^ vr (2) 

Proposition 2 A non barotropic and non isoenergetic divergence- free per- 
fect fluid energy tensor verifying ^ represents the local thermal equilibrium 
evolution of the ideal gas with specific internal energy e, temperature T , 
matter density r, and specific entropy s given by 

e(p,p) =e(^) = e(7r) - 1, T{p,p) = T{-k) = -e{7r) , (3) 

rip,p) = -r^, s{p,p) = kln , (4) 

e(7r) p 

e(7r) and /(vr) being, respectively, 

e(7r) =eoexp{ / V(vr)(i7r}, ^'(vr) = — — — -— ; (5) 

J (x(7r) -vr)(7r + l) 

/(vr) = /o exp{ / <t>{7r)d7T} , (/.(vr) = -r^ . (6) 

The above two propositions provide a complete algorithm, in four steps, 
to detect and characterize, in any given family of divergence-free perfect 
energy tensors T = {T = [u°'{x^), p{x^),p{x^)]}, those that represent an 
ideal gas evolving in local thermal equilibrium: 



step 1: Calculate the coordinate dependence of the space-time functions 

p/p = tt{x^) and u{p)/u{p) = xi^^) 

step 2: Determine the ideal gas subset of T by imposing the ideal gas 
hydrodynamic condition 0: dx A dvr = 0. 

step 3: Obtain, in this subset, the explicit expression of the indicatrix func- 
tion: x = x(7r). 

step 4: Calculate, from x(^)) the generating functions e{'ir) and /(vr) given 
in (0) and ©, and obtain from them the thermodynamic variables by 
using © and ©. 

3 THERMODYNAMIC STEPHANI UNIVERSES 

The conformally flat perfect energy tensor solutions to Einstein equations 
with nonzero expansion were first considered by Stephani 8^ and are usu- 
ally called Stephani universes. They were rediscovered by Barnes ,15., who 
obtained them as the conformally flat class of irrotational and shear-free 
perfect fluid space-times with nonzero expansion. 

In order to generalize the cosmological principle, Bona and Coll I^ looked 
for space-times admitting an iso-invariant synchronization, and without any 
hypothesis on the energy tensor, they likewise found the Stephani universes. 

An iso-invariant synchronization amounts the existence of an irrotational 
and shear-free observer u with non zero expansion such that the induced 
metrics on the orthogonal hypersurfaces are of constant curvature. Then, 
there exists an adapted coordinate system, u = {l/a)dt, such that the line 
element takes the form: 

ds^ = -a'^dt^ + n^6ijdx'dx^ (7) 

where 

a = RdRlnn , 17 = -. ^^^^ ^ (8) 

l + 2b{t)-r + lK{ty 

R{t), b{t) and K{t) being five arbitrary functions of time. 

The metric (jZj) is a solution to the perfect energy tensor of the form (fT]). 
with energy density and pressure given by 

P=j-,{R' + K-4b') , P=-P-§^ (9) 

Moreover, the expansion of u and the curvature of the spatial metric are 
homogeneous and respectively given by 

e{t) = ^^0, ^(t) = ^{K - 4b') (10) 



The conditions under which Stephani universes describe the evolution 
of a fluid in local thermal equilibrium, were directly obtained by Bona and 
Coll in |l9j, and later considered by Krasihski et al. UOj. The results that 
we need from [S] are here summarized in the following lemmas. 

Lemma 1 The necessary and sufficient condition for a Stephani universe 
to represent the evolution of a fluid in local thermal equilibrium, is that it 
admit a three-dimensional isometry group. 

A particular class of evolutions in local thermal equilibrium is that of 
barotropic evolutions; in our case, they are easily determined: 

Lemma 2 A Stephani universe represents the barotropic evolution of a fluid 
iff it admits a six- dimensional isometry group (Friedmann- Robertson- Walker 
space-time). 

The following lemma states, in thermodynamical terms, that Stephani 
universes do not admit isometry groups of dimension 4 or 5. 

Lemma 3 The dimension of the maximal isometry group of a Stephani 
universe that represents the non barotropic evolution of a fluid is three. 

This paper deals with Stephani universes describing non barotropic evo- 
lutions of fluids in local thermal equilibrium. In what follows, they will be 
called for short thermodynamic Stephani universes. 

A canonical form for them has also been obtained in P|. From it, and 
matching the two space- like coordinates r and z in a sole variable, say w, a 
straightforward calculation leads to: 

Proposition 3 The metric of the thermodynamic Stephani universes may 
be written 

ds^ = -a'^dt'^ + n'^6ijdx'dx^ (11) 

where 

w 
a = Rdn In L , fi = — L 

2z 
= ^(0 = 2z (12) 

1 + b{t)w ' "" = 1 + |r2 

R{t) and b{t) being two arbitrary functions of time. 

Its symmetry group is spherical, plane or pseudospherical depending on 
£ to be 1, or —1 and the Friedmann- Robertson- Walker limit occurs when 
b = constant. 

Furthermore, the energy density and pressure are given by 

P=^(i?^ + e-462) , P = -P-^^ (13) 

R"^ 3 a 



4 IDEAL GAS STEPHANI UNIVERSES 

Among the thermodynamic Stephani universes, what are those correspond- 
ing to ideal gases? 

In order to determine them, we shall study the restrictions that our above 
ideal gas condition imposes on the functions R{t) and b{t) that character- 
ize the general thermodynamic Stephani universes (|11|) dJ, following step 
by step the algorithm presented at the end of section 2. 

step 1: Coordinate dependence of n and x 

From the second equation in H13() . a direct calculation leads to: 

?- = tt{R,w) = --1 (14) 

p a 

where a is given in (|12jl and a is defined by 

a = a{R)^-^ , (15) 

Now we look for the indicatrix function x- From (|14() and ()15|) . we have 
u{p) _ dnp _ 1 



X = 



dnip-^) = vr - —Or ( - ) (16) 



u{p) dnp Orp 

but, from the definition of a in l\12\\ it follows 

.„ N l + (b-Rb')w 
» = »(R,nj)^ \^^^ , (17) 

where x' stands for dx/dR for any x = x{R). Then, from (|14j) and (|17)) it 
results: 



^■K + l J \b' R{tt + 1] 
and substituting in (fTC)) . we obtain: 

Lemma 4 For i/ie thermodynamic Stephani universes, the thermodynamic 
variable vr = p/p and the indicatrix function x = u{p)/u{p) take the expres- 
sions: 

,„, 11, , .a R i , , rio - , . 

X = x(^,^)=vr + ---(7r+l)[^ + - + (7r + l-a)^ (19) 

o o a Qj Qi 

where a is given by (|15l) and the prime indicates derivative with respect to 
the variable R. 



step 2: Ideal gas hydrodynamic condition: dx A dir = 

Expression (|19|) shows that, in agreement with the local thermal equilibrium 
condition, the indicatrix function x depends on the only energetic variables 
p and p, X = xiPiV)- Indeed, (|19jl can be written: 

X = X(vr,fi) = vr + ^ + ^(^ + 1)[(^ + l)Ai{R) + A2{R)] (20) 



where vr = p/ p and where 



Rh" A (f>\- ^^" «'^ 1 



MR)^--^, MR)^-jT--r-- (21) 



Thus, p being an effective function of R, the functions Ai and A2 can be 
considered as depending on p. Every choice of these two functions deter- 
mines different indicatrix functions that correspond to different media. How 
we must take Ai and A2 in the ideal gas case? 

The thermodynamic Stephani universes describing non barotropic evolu- 
tions of fluids in local thermal equilibrium are such that the function it{R, w) 
in H18() satisfies S^vr 7^ 0. Consequently condition dx/\dTT = is equivalent 
to dfixi'^: R) = 0. From (|2U|) this equation may be written: 

A[{R){tt + 1) + A'^iR) = (22) 

with Ai{R) and A2{R) given by (EIJ). 

But, TT and R being independent variables, from (|22|) it follows that the 
functions in (|2H) are constant, Ai{R) = Ci, A2{R) = C2. From now on, Ci 
and C2 will be called the principal constants. 

Moreover we can eliminate the function b{R) in the second of equations 
(|21() by using its values from the first one. Thus, in a first step, we obtain 
two equations for the two functions b{R), and p{R), this last one being 
related to a{R) by ((T3|) . In a second step, we put every solution b{R), p{R) 
for these equations in the expression (|13j) of the energy density p in order 
to obtain an equation for the expansion factor i?(t), which determines then 
the corresponding Stephani universe. More precisely, we can state: 

Theorem 1 A thermodynamic Stephani universe ()11() H12|) represents an 
ideal gas if, and only if, the metric function b{R) and the energy density 
p{R) satisfy the two equations: 

Ra' = -a(cia2 + c^a + 1) , b"R = -c,aH' (23) 

where a(R) is given in (fT5|l and the principal constants Ci and C2 are arbi- 
trary. 



Then, in terms of b{R) and p{R), the expansion factor R(t) satisfies the 
generalized Friedmann equation: 

p{R) = 1-^[R^ + e - Ah\R)] (24) 

A first glance at equations (|23|) (PU allows us to know that the ideal gas 
Stephani universes are, generically, a family of solutions depending on seven 
parameters, namely: the two principal constants Ci and C2 (we will see below 
that they determine five classes with different thermodynamical properties), 
the initial values Oq, Po and Rq for the functions a{R) p{R) and R{t) and, 
finally, two integration constants 61 and 62 giving the general solution h[R) 
to the second equation in (P5|l. In next section we study in detail equations 
(|23j) (|2l|) and we will analyze the properties of the models that we obtain 
depending on these seven parameters. 

step 3: Indicatrix function: x = x(^) 

Now we look for the explicit expression of the indicatrix function. Assuming 
equations ((231), the indicatrix (jT^ becomes a second degree polynomial in 
TT whose coefficients depend on the principal constants Ci and C2. More 
precisely, we have: 

Proposition 4 The indicatrix function x(^) of the ideal gas Stephani uni- 
verses is of the form: 

x(^) = /3vr2+77r + <5, (25) 

where P, 7 and 5 are constants related to Ci and C2 hy 

/3 = |, 7 = l + ^(c2 + 2c0, <5 = ^(c,+C2 + l). (26) 

step 4: Thermodynamic variables 

Finally, the expression (|25)) of the indicatrix function xi^^) determines, ac- 
cording to proposition 2, the other thermodynamic variables, and so, the 
thermodynamical properties of the fluid. Expressions Q and (jlj) show that 
the thermodynamic variables are determined by the generating functions 
e(7r) and /(vr) given by © and @ respectively. Consequently, the thermo- 
dynamic scheme depends on the reduced indicatrix x'- 

x(vr) = x(vr)-7r = /37r2 + 77r + <5, (27) 

with 7 = 7 — 1, and /?, 7 and 5 given by ()26() . Then, depending on the 
degree and the roots of x('7r) we find five different classes of ideal gas Stephani 
universes and, for every class, the integrals © and © admit simple analytic 
expressions. We summarize the results in the following: 

10 



Proposition 5 From a thermodynamical point of view and according to the 
values of the principal constants c^ and C2, the ideal gas Stephani universes 
may belong to five classes. For every one of them, the reduced indicatrix 
x(vr) and the generating functions e(7r) and /(vr) are given by 

Class 1 (ci = Ca = 0): xC^r) = 5 

/(tt) = /o expjSvr} , e(7r) = - — -—^ exp{37r} (28) 

Class 2 {c, = 0, c^ / 0): x(vr) = Mic^iir + 1) + 1] 



/(7r) = /o[(c.(7r + l) + l] 



e(vr) = ^-^[(c.(vr + l) + lf(^+i) 

Class 3 (A = c^ - 4c, =0, Ci / 0): x(7r) = ft [(c2(vr + 1) + 2]^ 

-12 



(29) 



/(tt) = /oexp 



C2[c2(7r + l)+2] 



, , [c,(7r + l) + 2]3 f -6(c2 + 2) 

e(7r) = eo 7 -to exp ■ 



(30) 



(7r + l)3 ^lc2[c2(7r + l) + 2]. 
Class 4 (A = c^ - 4ci > 0, c,^0): xl^r) = f (vr - 7r+)(7r - 7r_) 



3 



/(-) = /o (^^^^) ^ , vr,.-l + -L(-c,±VA) 
(^-7r+)^+(,,_^_)A- _3r c, + 21 

"^"^ = '° ^TTp ' ^^ = 2['^^^. 

Class 5 (A = c^ - 4ci < 0): ^(vr) = ^[civr^ + (2ci + C2)7r + c, + Ca + 1] 

„, X , f 2 2ci(7r + l) + ca \ 
; TT = fo exp < — ; arfftan ; > 

. ^ [X(vr)]^/^ f Ca + 2 2c,(7r+l)+ca ) 

The complete space-time ideal models, which will be obtained in next sec- 
tion by integrating equations (|23|) . depend on seven parameters. The above 
proposition shows that, among these seven parameters, the principal con- 
stants ci and C2 are the two thermodynamical ones. 
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5 IDEAL GAS STEPHANI MODELS 

The above section has obtained and explored the indicatrix function x foi^ 
ideal gases in Stephani universes, and constructed the generating functions 
/(vr) and e(7r) that directly give the thermodynamic variables. Here we 
complete the study by the exploration of the conditions (|23() for the existence 
of the indicatrix, which involve the metric functions. 

The study of these conditions (|23j) leads to different models that depend, 
generically, on seven parameters as pointed out above. The principal con- 
stants Ci and C2 distinguish the five classes of ideal gas Stephani universes 
already obtained in proposition |SJ But the analysis of the equations also 
distinguishes two families of models defined by the fact that the function 
a{R) given in ()15|) may be constant or not. Now we consider these two cases 
separately. 

5.1 Singular models: a'{R) = 

If a{R) = a,) 7^ 0, the first equation in H23() implies that the parameter 
ay of the model depends on the thermodynamic parameters, the principal 
constants Ci and C2. In these cases, we can easily integrate (fTC)) with a{R) = 
Uq to obtain p{R). On the other hand, the second equation in (|23j) leads to 
b' = Ci?~^^"Oj C being a constant. Then we obtain: 

Theorem 2 In the singular ideal gas Stephani models, a{R) = Oq, 
i) the energy density is given by: 

/ E> \ 3ao 

piR) = Po (^ j (33) 

where po and Rq are constants, and 

a) the metric function b{R) is given by one of the following expressions 

b{R)=b, + b,(^-^y if /i/O (34) 

b{R)=b, + bAn(^) if /i = . (35) 

where fi = 1 — Cia^, and bi, 62 7^ are arbitrary constants. 

As a consequence of (|23j) . the relationship between the constants Qqi P 
and Ci, C2 is given by: 

ci = 4(l-^) ' C2 = -(/i-2) (36) 



12 



These expressions show that both, a,, and ^u, are thermodynamical param- 
eters. Moreover, one can easily see how the above cases are related to the 
five classes of ideal gas Stephani universes established in proposition [51 

When yU = one has A = Cj — 4ci = 0, and these cases belong to class 
3. When ^ = 1 one has Ci = 0, Cj / 0, and they belong to class 2. All the 
other singular cases belong to class 4, because A = (/i/oo)^ > 0. Thus, if 
we take into account proposition El and (|H6|) we can state: 

Proposition 6 The singular ideal gas models, a{R) = Og, given in theorem 
121 belong to the following classes: 

if fi = 1, to class 2, with c-^ = . 

if fi = 0, to class 3, with c^ = — — ■ 

if ^ fi ^ I, to class 4, with Ci and C2 given by (|5^ . 

5.2 Regular models: a'{R) 7^ 

Equations ()23() are equivalent to equations (|21|) with Ai = c^ and A2 = C2] 

differentiating the first of these equations and considering the second one, 

one finds: 

a" 2a' R b" 

= 37 

a' a b' 

From this equation one can obtain 6 as a function of a. But this relation 

and the first equation in (|23|) imply the second one. Thus, we have: 



Proposition 7 The ideal gas Stephani universes with a'{R) 7^ are those 
for which 

Ra' = -afcio^ + C2a + 1) , b(R) = — ^ + b, (38) 

a{R) 

where Ci, C2, 61 and 62 7^ are arbitrary constants. 

Consequently, for regular ideal gas models, one needs to solve the first 
order differential equation (|38|) for the function a{R). The solutions can 
be obtained in an implicit form h(a, R) = depending on the degree and 
roots of the polynomial q{a) = CiO^ + C2a + 1. But an implicit expression 
for a{R) does not allow to determine explicit expressions for p{R) (solution 
to (|T5|) ) and b{R), and these expressions are necessary in order to establish 
the Friedmann equation (|24() and to obtain the metric tensor by solving it. 
We can overcame this shortcoming by a change of variables. Indeed, a been 
an effective function of R, the equation (jlSf) can be written 

1 dp 3a dR 3 

pda R da c^a^ + 020 + 1 ' 

13 



and consequently, the energy density takes the expression 

p{a) = p,Q\a) (40) 

where 

Q{a) = exp \ j —— \ , q{a) = c^c? + Cja + 1 (41) 

On the other hand, the first equation in (jSHJ takes the form aq{(i) dR/da = 

R, and its solutions depend on a simple function of q{a) and Q{a). Finally, 

Friedmann equation (PU may be considered as a differential equation on 

a{t) by changing: 

• dR R / . N 

R= a = ^a (42) 

da aq[a) 

All these considerations lead to the following: 

Theorem 3 The regular ideal gas Stephani models, a'{R) ^ 0, are those 
for which the metric functions R and b are of the form 

R{a) = ^ J\q{a)Q-^{a)\ , b{a) = ^ + b, (43) 

a * a 

where Q{a) and q{a) are given by 

Q{a) = exp I / -A- [ , q{a) = c^a^ + Caa + 1 (44) 

and a = a{t) is a solution to the generalized Friedmann equation 

c-i, C2, 61 and 62 7^ being arbitrary constants. 

The analytic expression of the function Q{a) depends on the polynomial 
q{a). Thus we can distinguish five different models depending on the values 
of the constants Ci and C2 . These cases coincide with the five classes of ideal 
gas Stephani universes established in proposition|5J The integration of these 
five cases leads to simple analytic expressions for the function Q{a): 

Proposition 8 According to the class to which the ideal gas Stephani uni- 
verse belongs, the functions q{a) and Q{a), that determine by (j43j) the metric 
function R when a'{R) 7^ 0, are given by: 

Class 1 (ci = Ca = 0): q{a) = 1 

Q(a) = exp{a} (46) 
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Class 2 


(Cl = 


^0, 


c.^ 


0): 


g(a) = ( 
Q{a) = 


"20.+ 1 

(c^a + l)- 


L/C2 




Class 3 


(A^ 


-cl 


-4ci 


= 0, 


Ci/0): 


q{a) = 


= i(c2a 


+ 2)2 










Q{a) = exf 


J ^ 


1 






IczCcza 


+ 2) J 




Class 4 


(A = 


-c\ 


-4ci 


>o, 


Ci^O): 


q{a) = 


= Ci(a- 


- aj^){a — a_) 






Q{a) = 


/a- 


a+ ^ VA 


, a± = 


i<- 


=2 ± ^)] 




\a- 


a J 


Class 5 


(A = 


- ^2 


-Ac, 


<0) 


: q{a) 


= CiO? + Caffl + 1 





(47) 



(48) 



(49) 



^/ N f 2 3(2cia + C2)l 

(^(aj = exp < —j=^ arg tan j=^ — > . (50) 

6 APPROXIMATE CLASSICAL IDEAL GAS 

Until now we have considered generic ideal gases characterized by the equa- 
tion of state p = krT. In this section we consider classical ideal gases which 
also satisfy the energetic equation of state e = c^T, and we analyze when 
a Stephani universe represents a classical ideal gas in local thermal equilib- 
rium. This energetic equation of state restricts the indicatrix function xi"^) 
of the ideal gases. More precisely, we have shown elsewhere jTj: 

Proposition 9 The indicatrix function of a classical ideal gas takes the 
expression: 

X(-) = i^ (51) 

where 7a = 1 + — ; is the adiabatic index, 1 < 7^ < 2. 

If a perfect fluid has the indicatrix function (|51() , then it represents a classical 

ideal gas with specific energy: 

e(^) = ^°:^ (52) 

7a - 1 - TT 

Prom this result and proposition |1] we can state: a thermodynamic 
Stephani universe never represents in a exact way a classical ideal gas in 
local thermal equilibrium. Nevertheless we should take into account that, in 
a relativistic framework, the classical ideal gas approximation works for low 
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temperatures. Thus, we can look for ideal gas schemes that satisfy the en- 
ergetic equation of state of a classical ideal gas at first order in T. This fact 
can be analyzed by studying the expressions ^ and their first and second 
derivatives for small T. Then, taking into account ((5J and proposition 9, we 
obtain 

Proposition 10 The necessary and sufficient condition for an ideal gas to 
have, at first order in the temperature, the energetic equation of state of a 
classical ideal gas, e(T) = CyT + o{T'^), is that the indicatrix function xi'^) 
and the specific energy e(7r) satisfy: 

X(0)=0, x'(0)=7a; e(0) = l (53) 

Now we can study the compatibility between (|53|) and the Stephani ideal 
gas indicatrix function (|25() at first order in vr, and we easily obtain 6 = 
and J = 'ja- Moreover, this means that the adiabatic index 7 fixes the 
parameters Ci and C2 of the classical ideal gas models. More precisely we 
have: 

Proposition 11 If an ideal gas Stephani universe represents a classical one 
at first order in the temperature, then the indicatrix function can be written: 

2 

x(^) = (7--)7r2+77r, 1<7<2 (54) 

and the principal constants are the following functions of 'j: 

Ci = 37 - 2 > 1 , C2 = 1 - 37 < -2 (55) 

This proposition implies that, if we impose on the thermodynamic pa- 
rameters ci and C2 the conditions ()55() . the ideal gas indicatrix function 
X = x(7r) satisfies (|SS)l. Then, only the models of class 4 in proposition |S1 
are admitted in this case. Indeed, from (|55|) we obtain: 

Ci/0, A = (ci-l)2 = 9(7-l)2 >0 (56) 

which agree with the conditions of proposition El that define the ideal gas 
Stephani universes of class 4. Moreover, if we determine the parameters 7r± 
and X± of the scheme (|31|) in terms of the adiabatic index 7 and we fix the 
constant cq in order to have e(0) = 1, we can state: 

Proposition 12 For the Stephani universes, the ideal gas thermodynamic 
schemes which approximate, at first order in the temperature, a classical one 
with adiabatic index 7, are generated by the functions 

1 
/(vr) = /o 



vr 



(7 - 2/3)7r + (7 - 1) 



7-1 

, efvr) 



1 + 



3(7-l)(7r + l) 



(57) 
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6.1 Classical singular models 

As a consequence of proposition IHl the singular models of class 4 satisfy 
7^ /I 7^ 1. Then, the relations (|36() and the restrictions ()55() lead to either 
Oo = 1 and fj, = 3(1 — 7) < 0, or a^ = 1/(87 ~ 2) and fi = 1 — Uq > 0. Then, 
if we take into account theorem |21 we can state: 

Proposition 13 There are two families of singular models that approxi- 
mate, at first order in the temperature, a classical ideal gas with adiabatic 
index 7; 

Models a,, = 1; 

Ki?)=A,(^) , HR) = W + h. i^^j (58) 

Models oq / 1/ 

/ Rn \ 7^273 / i? \ 7^^ 

p(R)=Po[^) ' ^(^) = ^^ + ^^ (£) (^^^ 

6.2 Classical regular models 

Proposition EJ and relations (|56() imply that the models of class 4 are the 
sole regular models compatible with the classical ideal gas approximation in 
question. If we write the parameters of the model in terms of 7 by using 
(|55j) . we obtain the following result: 

Proposition 14 For the regular models in theorem IS] that approximate, at 
first order in the temperature, a classical ideal gas with adiabatic index 7, 
the functions q (a) , Q{a), and R{a) take the expression: 



qia) = (a - 1) [(87 - 2)a - 1] , Q(o) 



R{a) = ^ 
a 



[(37 - 2)a - 1] 



37-2 



a-1 




3(7-1) 

7 




(87 - 2)a - 

1 

3(7-1) 


-1. 


(60) 



a-1 

7 ON THE PHYSICAL BEHAVIOR OF THE MODELS 

In this work we have presented Stephani models that describe evolutions 
in local thermal equilibrium of a generic ideal gas (section 5) or a classical 
ideal gas (section 6) . The space-time domain where these models have a good 
physical behavior are defined by complementary physical requirements. 
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Firstly, we have the Plebahski energy conditions that, for a perfect fluid, 
state —p < p < p ■ In terms of the hydrodynamic variable vr = p/p these 
conditions take the form — 1 < vr < 1 . Nevertheless, in the ideal gas case it 
is reasonable to consider positive pressures p. Thus we should impose on vr: 

< vr < 1 (61) 

On the other hand, the relativistic compressibility conditions J16j should 
be also required. Elsewhere ^SI we have shown that these conditions can 
be written in terms of hydrodynamic variables by means of the indicatrix 
function xip^v)-, aiid that for an ideal gas they become [T^ : 

^ / / N / 1 /^ N ^ 2x(vr)(l-x(7r)) 

<x(7r)<l, X[.'^)> --r-r\ T7 — T-jT (62) 



2vr + l ^' ' ^ ' ' (x(vr)-vr)(^ + l) 

These conditions imposed on the indicatrix function ()25() of the ideal gas 
Stephani universes imply the restriction of the parameters Ci and C2 which 
determine the thermodynamical properties of the models. Moreover, given 
a pair of values (ci, C2) physically compatible, the compressibility conditions 
H62() only hold for values of vr on a subset of the interval ]0, 1]. This means 
that the corresponding thermodynamic scheme is only well defined for these 
values of vr. Thus, for example, for a model of class 1 (that is, Ci = C2 = 0), 
conditions ^^ hold if vr G ]0, |[. 

For a classical ideal gas the compressibility conditions (|62() hold for every 
vr G ]0, 1[ if we assume that the adiabatic index 7 G ]1,2[. Nevertheless, 
for the Stephani models studied in section 6 which approximate a classical 
ideal gas at first order in the temperature, these conditions only hold for 
vr G ]0, vrm [, where the maximum value of vr depends on the adiabatic index as 
vrm = , ^ < 1 . It is worth remembering that this approximation 

7+V7' +4(7-2/3) ^ ^^ 

is valid only for small values of vr. Thus for vr near of the maximum value vr^, 
these models have a good physical behavior but they do not approximate a 
classical ideal gas. 

Once obtained the admissible values of the hydrodynamic variable vr for 
the different classes, we can look for the space-time domains where every 
model is physically reasonable. We should analyze the expression (|18|) that 
gives vr in terms of the spatial coordinate to and the metric function R{t). It 
is worth remarking that cj, as defined in p2j) . is a bounded coordinate in the 
case of spherical symmetry (e = +1) and, otherwise, it is not bounded. On 
the other hand, depending on the behavior of R{t) we could have models only 
valid at early or at present times. But all these analysis will be considered 
elsewhere. 

Finally, we want to comment about the relationship between the results 
by Sussman in ^J and ours. His models approximate a classical ideal gas 
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in some specific cases that, in our notation, correspond to the sphericaUy 
symmetric (e = 1) singular case (|58() with Cq = 1 and 61 = 0. In spite 
of this restriction, Sussman's paper presents a wide analysis on the space- 
time domains and observational features concerning the cases of a classical 
mono-atomic ideal gas (7 = 5/3) and a matter-radiation mixture (7 = 4/3). 
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